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Infrared divergences of QCD scattering amplitudes can be derived from an anomalous dimension matrix, which 
is also an essential ingredient for the resummation of large logarithms due to soft gluon emissions. We report 
a recent analytical calculation of the anomalous dimension matrix with both massless and massive partons at 
two-loop level, which describes the two-loop infrared singularities of any scattering amplitudes with an arbitrary 
number of massless and massive partons, and also enables soft gluon resummation at next-to-next-to-leading- 
logarithmic order. As an application, we calculate the infrared poles in the qq — >■ tt and gg — >■ tt scattering 
amplitudes at two-loop order. 



1. Introduction 

Due to the complicated self-interactions among 
the gauge bosons, the structure of infrared (IR) 
singularities is a highly non-trivial property of 
non-abelian gauge theories such as quantum chro- 
modynamics (QCD). The knowledge of IR diver- 
gences in the scattering amplitudes is essential 
to the derivation of factorization theorems [1]. 
And while factorization theorems guarantee the 
absence of IR divergences in sufficiently inclusive 
observables, in many cases large Sudakov loga- 
rithms remain after this cancellation. A detailed 
control over the structure of IR poles in the vir- 
tual corrections to scattering amplitudes is a pre- 
requisite for the resummation of these logarithms 
beyond the leading order ^2, 3 • 

For scattering amplitudes involving only mass- 
less partons, Catani was the first to predict the 
singularities at two-loop order apart from the sin- 
gle pole [4], whose general form was only un- 
derstood much later in [5H2- In recent works 
[8l|9], it was shown that the IR singularities of 
on-shell amplitudes in massless QCD can be de- 
rived from the ultraviolet (UV) poles of opera- 
tor matrix elements in soft-coUinear effective the- 
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ory (SCET) [TMl]- They can be subtracted 
by means of a multiplicative renormalization fac- 
tor, whose structure is constrained by the renor- 
malization group. The corresponding anomalous- 
dimension matrix is constrained by soft-collinear 
factorization, the non-abelian exponentiation the- 
orem, and the behavior of scattering amplitudes 
in two-parton collinear limits. With these con- 
straints, the results of [JHZ] can be easily under- 
stood, and it was further proposed that the sim- 
plicity of the anomalous-dimension matrix holds 
not only at one- and two-loop order, but may 
in fact be an exact result of perturbation theory. 
This possibility was raised independently in [13], 
while potential new structures at three-loop were 
investigated in [M] . 

It is relevant for many physical applications 
to generalize these results to the case of massive 
partons. The IR singularities of one-loop ampli- 
tudes containing massive partons were obtained 
some time ago in |15| . but until very recently 
little was known about higher-loop results. In 
the effective theory language used in [Hlin], one 
needs heavy-quark effective theory (HQET) (see, 
e.g., [16] and references therein) besides SCET, 
and the simplicity of the anomalous dimension 
matrix observed in the massless case no longer 
persists in the presence of massive partons. While 
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the non-abelian exponentiation theorem stih re- 
stricts the allowed color structures, important 
constraints from soft-coUinear factorization and 
two-parton coUinear limits are lost. In |17j , the 
generic form of the anomalous-dimension matrix 
are given, with two universal functions -Fi and 
/2 representing three-parton correlations left un- 
specified. 

Very recently, in [TSKH], the two functions Fi 
and /2 were calculated analytically, and there- 
fore the generic two-loop anomalous-dimension 
matrix with arbitrary number of massive and 
massless partons is now completely determined. 
These results provide a general, unified, and an- 
alytic description of the IR singularities arising 
up to two-loop order in arbitrary scattering am- 
plitudes with massive and massless external legs 
in unbroken gauge theories such as QCD and 
QED. They also form the basis for the systematic 
resummation at next-to-next-to-leading logarith- 
mic (NNLL) order of large Sudakov logarithms 
affecting the rates for such processes. As a first 
application, the IR poles in the two-loop virtual 
amplitudes for qq tt and gg — > tt processes 
were determined in [THUH]. Utilizing these re- 
sults, an approximate next-to-next-to-leading or- 
der (NNLO) formula and an NNLL resummation 
formula for top quark pair production at hadron 
colliders were later derived in [20l[2T] with de- 
tailed discussions on numerical results and phe- 
nomenological implications. In this talk, we first 
briefly review the generic setup in [8l|9l[T7] and 
then present the results of [l8llT9] . 



2. Infrared singularities and anomalous di- 
mension matrix 

We denote by | A^n(e, {p}, {m})) a UV- 
renormalized on-shell rt-parton scattering ampli- 
tude with IR singularities regularized in d = 4— 2e 
dimension. Here {p} — {pi, . . . ,p„} and {m} = 
{rni, . . . , m„} denote the momenta and masses of 
the external partons. The amplitude is a func- 



tion of the Lorentz invariants s.; 



and pf — 



2aijp.rpj +iO 



where the sign factor cTy = +1 



if the momenta Pi and pj are both incoming or 
outgoing, and aij = — 1 otherwise. For massive 
partons we define four velocities w,; — Pi/rrii with 



vf = 1. We further define the recoil variables 
Wij = —<7ijVi ■ Vj — iO. We use the color-space 
formalism in [22] extended to the case with Wil- 
son coefficients and effective operators, which has 
been explained in detail in ^19^21^ and will not be 
repeated here. 

It was shown in [8j[9l[T7] that the IR poles of 
such amplitudes can be removed by a multiplica- 
tive renormalization factor Z^^(e, ijpj, |m), u), 
which is the same one used to renormalize the 
UV divergences of the operator matrix elements 
in the effective theory. The renormalization fac- 
tor satisfied the renormalization group equation 
(RGB) 

Z (e, {p}, {m}, fi) = -r({p}, {m}, p) . 



c?ln/x 



(1) 



In the case of only massless partons, the anoma- 
lous dimension matrix F has the simple form [8ll9 



where the absence of three-parton correlation 
terms was explicitly demonstrated in [6l[7]. In 
[S1[S1[T3], it was conjectured that this result may 
hold to all orders of perturbation theory. On the 
other hand, when massive partons are involved in 
the scattering process, then starting at two-loop 
order, correlations involving more than two par- 
tons appear. At two-loop order, the general struc- 
ture of the anomalous dimension matrix is |17) 
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Tcusp {^^S 



E ^7cusp(/3/j,as) + E'^^("«) 
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-Slj 



r'TfT'jT-tcFi{Pij,PjK,PKi) 



(I,J.K) 

+ E E T}T'}n /2 ( /3/.7, In ■ 



-crjk vj ■ Pk 
~<Jik VI ■ Pk , 

(2) 



where the lower case indices i, j, k denote mass- 
less partons and the upper case indices /, J, K de- 
note massive partons in the external states. The 
cusp angles /J/j are defined via cosh/?/,/ = wij. 

The coefficient functions 7cusp(as), Yi'^s) (for 
i = q, g) have been determined to three-loop or- 
der in '9^ by considering the case of the massless 
quark and gluon form factors 23: . Similarly, the 
coefficients 7^ (as) for massive quarks and color- 
octet partons such as gluinos have been extracted 
at two-loop order in (17j using the known two- 
loop anomalous dimension of heavy-light currents 
in SCET [23]. In addition, the velocity-dependent 
function 7cusp(/?j Q^s) has been derived from the 
known two-loop anomalous dimension of a cur- 
rent composed of two heavy quarks moving at dif- 
ferent velocity 25 , 26 . The most difficult ones to 
compute are the two functions Fi and /2, which 
start at two-loop order and involve correlations 
among three partons. They have been derived 
in closed analytic form in [T8l[T9], and were con- 
firmed in the non-physical region by a numerical 
result in p7l[28]. 

It is worth noting here that the anomalous di- 
mension matrix T not only determines the IR di- 
vergences in the amplitudes, but also enters the 
renormalization group evolution of the hard func- 
tions in the factorization formulas, and therefore 
serves as the basis for predicting and resumming 
logarithmic enhanced terms to all orders in as- 

3. The functions Fi and /2 

The calculation of Fi and /2 involves the eval- 
uation of the Feynman diagrams in Figure [T] in 
the soft limit. More formally, we evaluate the 
two-loop vacuum matrix element of the operator 
Os — Sv^ Sy^ Sy^ , which consists of three soft Wil- 
son lines along the directions of the velocities of 
the three partons, without imposing color conser- 
vation. For Fi, all the three directions are time- 
like; while for /2, one of them become light- like. 
We will describe the calculation of Fi and then 
obtain /2 from a limiting procedure. 

The operator Os is renormalized multiplica- 
tively so that OgZs is UV finite. The anoma- 
lous dimension of this operator, which equals 
Gir^'^T^T^T^ Fi{pi2, [323, can then be ob- 




Figure 1. Two- loop Feynman graphs (top row) 
and one-loop counterterm diagrams (bottom row) 
contributing to the two-loop coefficient of the 
renormalization factor Z^. 



tained from Zg by using Eq. ([T]). To determine 
Zg, we need to separate the UV and IR diver- 
gences in the amplitude, for that we introduce 
an offshellness for the heavy quarks to regulate 
the IR divergences, while the UV divergences are 
regulated dimensionally. 

The planar and counter-term diagrams in Fig- 
ure [Dean be evaluated with standard techniques. 
The contribution to Fi from those diagrams reads 



(2) planar+CT 



^1^ 



X coth /3_ 



eijK (3ki coth Pki 



I,J,K 



ij 



/3/2/ + 2/3/jln(l-e-2^") 
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The diagram involving the triple-gluon vertex is 
the most technically challenging one. We have 
computed this diagram using a Mellin-Barnes 
representation with the result [TH] 



^(2)„o„-pla„ar^_4 ^ ,^^^^2^^^^e0th 
I,J,K 



For Euclidean velocities, our result agrees numer- 
ically with a position-space-based integral repre- 
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sentation derived in [2T. Combining all contri- 
butions, we find 



I,J,K 



(3) 



where 

r(/3) = /3coth/3, 



g{(3) = coth^ 



+ 2^1n(l 



Li2(e-^^) + y 



(4) 



The function /2 can be obtained from the above 
result by writing W23 = -(T23V2 • Pz/m^-, w^i = 
—a^iVi ■ Pz/tti^ and taking the limit 7713 — at 
fixed Vi ■ P3. In that way, we obtain 



-Cr23V2 ■ P3 



..ONI -0^23^'2 -^3 f^s 

= -4g(A2) In . (5) 

Whether a factorization of the three-parton terms 
into two functions depending on only a single cusp 
angle persists at higher orders in Ug is an open 
question, but our results certainly suggest a hid- 
den structure, which seems worthy of further ex- 
ploration. 

It is interesting to investigate the behavior of 
the two functions Fi and /2 in the limit where two 
velocities coincide. This is relevant, e.g., for the 
situation when two heavy quarks are produced 
nearly at rest. Contrary to naive expectations 
based on antisymmetry, the two functions don't 
vanish in this limit. Instead, they tend to a finite 
value due to the presence of Coulomb singulari- 
ties. On the other hand, in the limit when all the 
parton masses are small, both functions vanish 
like {mimj)'^ /sij, in accordance with a factor- 
ization theorem proposed in p9ll30]. 



4. Top quark pair production 

Top quark pair production at hadron colliders 
is an important process in the standard model. It 



is also the simplest example which involves two 
heavy quarks where Fi and /2 can contribute. 
Therefore, as a first application, it is interest- 
ing to use our results in this case. In [151 [T5] 
we have worked out the relevant anomalous di- 
mension matrices for both qq — > ti and gg — tt 
partonic processes, and derived the IR poles in 
the two-loop amplitudes. For the qq channel, our 
results agree with the numerical ones from |31) . 
with the analytic results for some of the color co- 
efficients given in j32, 33j, and with the results in 
the small- mass limit from [S^. For the gg chan- 
nel, results in the literature are available only in 
the small-mass limit [35j , and we have checked the 
agreement of our exact results with this limiting 
case. Later at the RADCOR 2009 conference, 
M. Czakon claimed that they had confirmed our 
results in the gg channel [36]. 

Given the knowledge of the anomalous dimen- 
sion matrices, it is possible to derive approximate 
formulas beyond next-to-leading order for vari- 
ous differential cross sections as well as the total 
cross section for this process. This was done for 
the invariant mass distribution in [18] where all 
the threshold enhanced terms at NNLO were ex- 
actly predicted. In [T^], these threshold enhanced 
terms were resummed to all orders in as at NNLL 
accuracy, where also the total cross section was 
computed by integrating over the invariant mass. 
Details about these two works are discussed by 
A. Ferroglia in this proceeding. 

5. Conclusions 

The IR divergences of scattering amplitudes in 
non-abelian gauge theories can be absorbed into a 
multiplicative renormalization factor, whose form 
is determined by an anomalous-dimension matrix 
in color space. For processes with only massless 
partons in the external states, a simple form of 
the anomalous-dimension matrix has been con- 
jectured. Although a lot of supporting argu- 
ments were given, a rigorous proof is still lacking. 
When massive partons are involved, at two-loop 
order, the anomalous-dimension matrix contains 
pieces related to color and momentum correla- 
tions among three partons as long as at least two 
of them are massive. This information is encoded 
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in two universal functions: Fi, describing cor- 
relations among three massive partons, and /2, 
describing correlations among two massive and 
one massless parton. These two functions have 
been calculated analytically in [T8lfT9] , and there- 
fore the IR divergences of any two-loop scattering 
amplitude with arbitrary number of massive and 
massless external particles are completely under- 
stood. These results also provide the basis for 
a systematic resummation of Sudakov logarithms 
at NNLL order, and are relevant to a large class 
of interesting hadron collider processes and thus 
have important implications for precision mea- 
surements at the LHC. As first applications, we 
worked out the IR poles in the two-loop ampli- 
tudes for qq — ^ ti and gg tt processes, derived 
an approximate NNLO formula for the invariant 
mass distribution as well as total cross section 
for top quark pair production at hadron collid- 
ers, and also resummed the threshold enhanced 
terms to all orders in as at NNLL accuracy. 
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